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Mass and Spin of GRS 1915+105 Based on a Resonance Model of QPOs 
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Abstract 

We demonstrate that the four high-frequency QPOs observed in GRS 1915+105 can be interpreted as 
the oscillation modes on disks that non-lincarly and resonantly interact with a warp. The warp is assumed 
to be a low-frequency global pattern on the disk. This identification suggests that the relevant mass, M, 
and spin parameter, a, are M = 13-14M Q and a — 0-0.15. 
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1. Introduction 

GRS 1915+105 is the most extraordinary variable X- 
ray binary known. It exhibits a rich set of time variability 
in X-rays. So far, four high-frequency quasi-periodic oscil- 
lations (HFQPOs) have been observed in GRS 1915+105. 
Their typical frequencies are 41 Hz, 67 Hz, 113 Hz, and 
168 Hz (a review by McClintock, Remillard 2003). (These 
QPOs are hereafter called, in turn, the 40 Hz, 67 Hz, 113 
Hz, and 168 Hz QPOs, although the frequencies may not 
represent real ones.) The 113 Hz and 168 Hz QPOs are a 
pair in the sense that their frequency ratio is close to 2 : 3. 
GRS 1915+105 is one of three black hole binaries (the oth- 
ers are XTE J1550-564 and GRO J1655-40) with a pair 
of high-frequency QPOs whose frequency ratio is close to 
2 : 3. 

Time variabilities are good tools to investigate the 
structure of rotating gas in a strong gravitational field 
of black holes (e.g., Kato 2001). There are a few models 
proposed to explain HFQPOs. Abramowicz and Kluzniak 
(2001) and Kluzniak and Abramowicz (2001) suggest that 
the pair HFQPOs are due to a parametric resonance, and 
that their frequency ratio must be commensurable, e.g., 2 
: 3. The 67 Hz QPO, which is known to be rather stable, is 
suggested by Nowak et al. (1997) to be a trapped g-mode 
oscillation resulting from the presence of the maximum 
of the radial epicyclic frequency in a strong gravitational 
field (Okazaki et al. 1987; Perez et al. 1997). 

Recently, Kato (2004a, 2004b) proposed a resonance 
model of HFQPOs. He assumed that a warp, which is 
a low-frequency and global pattern, is present on ax- 
isymmetric disks. Some of the disk oscillations then 
nonlincarly and resonantly couple with the warp (Kato 
2004a, referred hereafter to Paper I), and are amplified 
or damped (Kato 2004b, referred hereafter to Paper II). 
The purpose of this letter is to examine, from the view- 
point of frequency, whether the observed four HFQPOs of 
GRS 1915+105 can be explained by this resonance model, 
and how much the resulting mass and spin of the central 
source actually are. 
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Fig. 1. Feedback processes of nonlinear resonant interac- 
tions acting on oscillations. The original oscillations are char- 
acterized by oj, m, and n. Since the warp corresponds to a 
wave mode of to = 0, m — 1, and n = 1, the nonlinear in- 
teraction between them brings about intermediate modes of 
oscillations of mil, and nil. To these intermediate 
oscillations, the disk resonantly responds at a certain radius. 
Then, the intermediate oscillations feedback to the original 
oscillations after the resonance. This feedback process ampli- 
fies or dampens the original oscillations, since resonances are 
involved in the interaction processes. (After Paper I) 
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2. Spin-Mass Relations Derived from QPO 
Frequencies 

We consider a disk that is deformed by the presence of 
a warp. The warp is a global one-armed (m = 1) pattern 
with a low frequency. When the disk is not deformed and 
axisymmetric, there are two kinds of oscillation modes, 
i.e., the p- modes and g- modes (e.g., Kato et al. 1998). 
Because of the deformation, the disk oscillations non- 
lincarly couple with the warp. Some of the coupled os- 
cillations (we call them intermediate oscillations) can res- 
onate with the disk at some particular radii, depending 
on the modes of the oscillations. These intermediate os- 
cillations can further interact non-linearly with the warp 
to feedback to the original oscillations (see figure 1). The 
feedback to the original oscillations passing through a res- 
onance is of interest, since by energy exchange with the 
disk through the resonance the amplitudes of the original 
oscillations change with time. 

Figure 1 demonstrates schematically the case where the 
warp has no precession. This case is considered first. 
Next, the case where the warp has a slow precession with 
a low frequency, u> p , is discussed. 

(a) The case where warp has no precession 

The analyses in Paper I show that there are two kinds 
of resonances, i.e., horizontal and vertical. The analyses 
further show that if we consider long-wavelength oscilla- 
tions, the resonances occur near to the radii where one of 
the following relations is realized: 

k = 0, V2Q± — Ok = Ok = 2k, V3Q,± — Ok = k,(1) 

where Ok(0 is the (relativistic) Kepler frequency and 
Oj_(r) and n(r) are the frequencies of the vertical and 
radial cpicyclic oscillations, respectively. The radii where 
one of the relations of equation (1) is satisfied are de- 
noted hereafter, in turn, as ro, r^/g, ?*4, and ryg, which 
are functions of the spin parameter, a. In the case of the 
Schwarzschild metric (a = 0), the radii satisfying the above 
relations are, in turn, 

ro = 3r g , ry2 = 3.62r g , r 4 = 4.0r g , = 6.46r g ,(2) 

where r g is the Schwarzschild radius, defined by r g = 
2GM/c 2 , M being the mass of the central black hole. 

The most interesting resonances are those at 7-4, since 
they are amplified (Paper 11). The resonances at are 
realized by horizontal resonances of g-mode oscillations 
(including the p-mode oscillations without node in the 
vertical direction), and the frequencies of the resonant os- 
cillations are ?ti*Ok ±K at r = n, depending on the modes 
of oscillations, where m* is zero or an integer (Paper I 
and II). In the case of one-armed (m = 1) oscillations, 
oscillations of 2(m*OK ± «) are also expected (Paper II). 
The most prominent oscillations among them arc those of 
2(Ok — k) and 20k — k (or Ok + k), which are 2k and 3k 
in the present case of no precession (Paper II). Hence, we 
think that the oscillations of 2k and 3k at r± are the ob- 
served 2 : 3 pairs of HFQPOs. That is, we take 3k at 



to be 168 Hz. This gives a relation between the spin pa- 
rameter, a, and the mass of the central object, M. which 
is shown in figure 2 by the curve labelled by 168 Hz. 

The resonant oscillations at other radii (k = 0, V20 j_ — 
Ok = and — Ok = k) are all damped when the 

warp is steady (Paper II). Among these damping oscil- 
lations, however, weakly damping ones will be amplified 
and observable when the warp is in a growing stage. In 
this sense, the resonant oscillations at the above radii are 
also worth investigating. Observations show occasional 
appearances of an ~ 40 Hz QPO and the so-called 67 Hz 
QPO. Considering that these frequencies are relatively low 
compared with the pairs, we assume that they are reso- 
nant oscillations at r /g. Then, they are all vertical reso- 
nances of g-mode oscillations (Paper I) , and their frequen- 
cies are 777*Ok±k at (Paper I), where m* is again zero, 
or a positive integer. Among the oscillations of m*OK ± k 
at r^, Ok — k is too low for our present interest. Hence, 
as the next lowest frequencies of resonant oscillations, we 
adopt k and 20k — «• That is, we regard the 40 Hz and 67 
Hz QPOs as resonant oscillations with frequencies k and 
20k — k at r = r^. Recent observations show that the 
frequency of the so-called 67 Hz QPOs is 69.2±0.15 Hz 
(Strohmayer 2001). Hence, we adopt here K r=r ^. = 40 Hz 
and (20k — K) r=r = 69 Hz. They give two independent 
relations between a and M, which are also shown in figure 

2 with labels of 40 Hz and 69 Hz. In order to demonstrate 
how much the curves shift on the diagram when the fre- 
quencies are changed, the curves of 41 Hz and 68 Hz are 
also shown in figure 2. 

Figure 2 shows that three curves of 168 Hz, 69 Hz, and 
40 Hz cross near to a point, although not at a point. As 
the curves of 68 Hz and 41 Hz show, curves move to the 
left of the diagram if the labelled frequency is increased. 
That is, if we adopt 39.9 Hz instead of 40 Hz, while keep- 
ing the other frequencies unchanged at 168 Hz and 69 
Hz, the three curves cross at a point, giving a ~ and 
M ~ 12.7M©. If we adopt 69.3 Hz instead of 69 Hz, while 
keeping the others unchanged at 168 Hz and 40 Hz, the 
three curves cross at an another point, giving a ~ 0.02 
and M ~ 12.9M Q . As another example, the frequency of 
the 67 Hz QPO is increased to 70 Hz, while keeping the 
others unchanged at 168 Hz and 40 Hz. Then, the rela- 
tions among the three curves become as shown in figure 
3. This shows that three curves cross at a point if the 
frequency 40 Hz is slightly increased to 40.2 Hz, while 
keeping the other frequencies unchanged at 168 Hz and 
70 Hz. This case gives a ~ 0.04 and M ~ 13.1M Q . Figure 

3 further shows that the three curves cross at a point if 
the frequency of the 168 Hz QPO is decreased to 167 Hz, 
giving a ~ 0.06 and M ~ 13.3M These considerations 
show that a unique set of a and M is possible for suitable 
sets of three frequencies near to (168 Hz, 70 Hz, and 40 
Hz). The mass obtained in this way is consistent with 
M = 14 ± 4Mq observationally derived by Greiner et al. 
(2001). 

(b) The case where warp has precession with frequency 
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So far, evaluations of the spin parameter, a, and the 
mass, M, are made under the assumption that the warp 
has no precession. When the central black hole has 
spin, however, the warp generally has a precession. If 
the warp is precessing with frequency u) p , the resonant 
radii are modified from those given by equation (1). The 
analyses in Paper I suggest that in the lowest order of 
approximations, the resonance conditions J1k = 2k and 
— Qk = k, for example, are changed, respectively, 

to 



2k ± uiry and V?>Q\ 



K i LOr 



(3) 



where both ± are possible, but we hereafter take the up- 
per sign (+), since this case has a better agreement with 
observations. The frequencies of the pairs at the radius 
of Qk = 2k + Up are then 217k — 2k (= 2k + 2w p ) and 
2J1k — k (= 3k + 2cl>p).^ They are denoted by u>\ and lo u , 
respectively The frequency difference is still k, but the 
ratio of the pair frequencies is deviated from 2 : 3. The 
frequencies corresponding to the 40 Hz and the 67 Hz os- 
cillations are now the frequencies of k and 2f2 — k at a 
radius of \/3± — = k + cj p . 



The case of 



0.5 Hz with 168 Hz, 69 Hz, and 39 



Hz is shown in figure 4. Since uj p = 1.5u>i — u) u , we have 
lo\ = 2u> u /3 + 2cj p /3. Hence, the adoption of oj u = 168 Hz 
and uj p = 0.5 Hz means that wc have adopted uj\ = 112.3 
Hz. In figure 4, the three curves cross near to a point. To 
make the three curves cross closer, keeping 168 Hz and 
69 Hz unchanged, the frequency 39 Hz should be slightly 
decreased as 38.9 Hz. Then, the crossing point gives a ~ 
0.11 and M ~ 13.8M Q . As an example of a larger w p , the 
case of u! p = 2.0 Hz is shown in figure 5 with 168 Hz (and 
175 Hz), 69 Hz, and 37.5 Hz. This figure shows that the 
curves of the 40 Hz QPO and of the 67 Hz QPO cannot 
touch each other unless the frequency of the 40Hz QPO is 
taken to be as low as 37.5 Hz when 70 Hz is adopted for the 
67 Hz QPO. (If a frequency lower than 70 Hz is adopted 
for the 67 Hz QPO, a frequency smaller than 37.5 Hz is 
required for the 40 Hz QPO.) The figure further shows 
that the frequency of the 168 Hz QPO must be higher 
than 175 Hz for the three curves to cross at a point. 

3. Discussion 



The results given in section 2 suggest that the precession 
of the warp must be low. That is, if lu p is as large as ui v = 
2.0 Hz, the requirement that the three curves cross at a 
point on the a-M diagram imposes frequencies somewhat 
different from observations on the 168 Hz, 67 Hz, and 40 
Hz QPOs. These considerations suggest that reasonable 



values of < 



and M are i 



0-0.7 Hz, a = 0-0.15, and 



M = 13-14M Q . 

Let us discuss the possibility as to whether trapped c- 
mode oscillations are the warps required here. Analyses 
by Silbcrgleit ct al. (2001) show that the frequency of the 
trapped c-mode oscillations is ~ 2 Hz when a ~ 0.1 and 



The case where the upper frequency is Ck + t(= 3k + oj p ) is also 
possible. 
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Fig. 2. Spin-mass relations derived from QPO frequencies, 
based on our resonance model. The warp is assumed to have 
no precession. Three curves crossing near to a point on the 
diagram are for 168 Hz, 69 Hz, and 40 Hz. In order to demon- 
strate how much the curves shift on the diagram when the 
frequencies are changed, the curves for 68 Hz and 41 Hz are 
also shown. If the frequency 40 Hz is slightly decreased to 
39.9 Hz, the three curves cross at a point, giving a ~ and 
M ~ 12.7-Mq. If the frequency 69 Hz is slightly increased to 
69.3 Hz, the three curves cross at an another point, giving 
a ~ 0.02 and M ~ 12.9M . 

M~ 14M Q , which seems to be slightly larger than the pre- 
cession frequency discussed in section 2. Furthermore, the 
trapped region of the c-mode oscillations is rather narrow 
in the radial extent (Silbcrgleit ct al. 2001). These re- 
sults suggest that the trapped c-mode oscillations are not 
directly related to the warp. GRS 1915+105 has many 
LFQPOs. Some of them might be related to the warp. 
The excitation of time-dependent low-frequency warps on 
disks is conceivable, if we remember that the global disk 
structure of GRS 1915+105 occasionally changes by a 
phase transition between hard and soft states. 

Nowak et al. (1997) suggested that the 67 Hz QPO is 
a trapped g-mode oscillation. In our present model, this 
is not the case, if we regard 113 Hz and 168 Hz QPOs 
as being 2 : 3 pairs whose frequencies are close to 17k 
and f^K + k at the radius r = r^. Adoption of this model 
implies that k at r — is ~56 Hz. When a = 0, this 
frequency is equal to K max (the maximum value of k), 
since is the radius where k = K max . This implies that 
when a = 0, the 67 Hz QPO cannot be trapped g-mode 
oscillations, since trapping requires that their frequencies 
must be smaller than K max - When a is large, the radius r& 
is outside of r max , and K max is larger than n(ri). However, 
calculations show that K max cannot become as large as 67 
Hz unless a ~ 1. This means that the model of the trapped 
g-mode oscillations and our present resonance model are 
incompatible. 
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Fig. 3. Spin-mass relations derived from the QPO frequen- 
cies, based on our resonance model. The warp is assumed to 
have no precession. The three curves are for f68 Hz, 70 Hz, 
and 40 Hz. These curves cross near to a point on the diagram. 
If the frequency, 40 Hz, is slightly increased to 40.2 Hz, the 
three curves cross at a point, giving a~0.04 and m~ 13.1A4q. 
Similarly, if the frequency of the 168Hz QPO is decreased to 
167 Hz, others being kept at 70 Hz and 40 Hz, the three curves 
cross at a point, giving a ~ 0.06 and M ~ 13.3Mq. 



Fig. 4. Spin— mass relations derived from the QPO frequen- 
cies, based on our resonance model. The warp is assumed to 
precess with a frequency of 0.5 Hz. Three curves are for 168 
Hz, 69 Hz, and 39 Hz. If the frequency of 39 Hz is slightly 
decreased to 38.9 Hz, the three curves cross closer, giving 
a ~ 0.11 and M ~ 13.8M©. 



In our resonance model, the frequencies of the resonant 
oscillations are not robust. This is because the resonant 
oscillations are not radially trapped, and can propagate 
in the radial direction, while changing the wavenumbers 
and frequencies. This will be one reason why the QPOs 
have large Q- values. Another reason for the large Q- values 
would be that different modes of oscillations with different 
coupling paths have close frequencies. If our resonance 
model represents real situations, some resonant oscilla- 
tions at y/2Q± — f^K = K may be expected to be observed. 

Finally, the azimuthal and vertical node numbers (to, n) 
characterizing the resonant oscillations are summarized. 
The pair oscillations at are p-mode or g-mode oscil- 
lations, coupling horizontally with the warp. The lower 
frequency oscillations of the pairs are the oscillations of 
m = 1 (one-armed oscillations). They are the p-mode os- 
cillation of n = or the g-mode oscillation with arbitrary 
n, as discussed in Papers I and II. The higher frequency 
oscillations of the pairs are (i) the g-mode oscillations of 
to = 1 or 2 with arbitrary n or (ii) the p-mode oscillation 
of to = 1 or 2 with n = 0. The 40 Hz QPO and the 67 Hz 
QPO at r = are due to the vertical resonance of the 
g-mode oscillations. The 40 Hz QPO is the oscillation of 
(to = 0, n = 2) and the 67 Hz QPO is the oscillation of 
(m = 2, n = 2). 

The author thanks Dr. Y. Kato for providing some 
references concerning observations. 
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Fig. 5. Spin— mass relations derived from the QPO frequen- 
cies, based on our resonance model. The warp is assumed to 
precess with a frequency of 2.0 Hz. The four curves are for 
175 Hz, 168 Hz, 69 Hz, and 37.5 Hz. 
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